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Abstract. The leading asymptotic large-scale behaviour of the spatially bipartite en¬ 
tanglement entropy (EE) of the free Eermi gas infinitely extended in multidimensional 
Euclidean space at zero absolute temperature, T = 0, is by now well understood. Here, 
we present and discuss the first rigorous results for the corresponding EE of thermal 
equilibrium states at T > 0. The leading large-scale term of this thermal EE turns out 
to be twice the first-order finite-size correction to the infinite-volume thermal entropy 
(density). Not surprisingly, this correction is just the thermal entropy on the inter¬ 
face of the bipartition. However, it is given by a rather complicated integral derived 
from a semiclassical trace formula for a certain operator on the underlying one-particle 
Hilbert space. But in the zero-temperature limit T ), 0, the leading large-scale term 
of the thermal EE considerably simplifies and displays a ln(l/T)-singularity which 
one may identify with the known logarithmic enhancement at T = 0 of the so-called 
area-law scaling. 
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1 . Introduction and main results. 

In recent years, entanglement entropy has turned out to be a useful and much 
studied quantifier of nonclassical correlations between subsystems in composite quantum 
systems [1]. In particular, given the (pure) ground state of a spatially large many- 
particle system and reducing (or localising) it to a spatial subregion fl, we denote the 
von Neumann entropy of the resulting (mixed) substate by S'(0, fl) and call it the local 
ground-state entropy. The spatially bipartite entanglement entropy (EE), defined for a 

t Published in slightly different form as J. Phys. A: Math. Theor. 49 (2016) 30LT04 (9pp). 

[There it should read (Tq/T)^/^ in Eq. (14b).] See also arXiv: 1501.03412v3 [quant-phj. 
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hounded as the (quantum) mutual information relative to the complement of and 
denoted by -ff(0, f2), then simply equals 25(0, ff) by the purity of the ground state. This 
(ground-state) EE quantihes, to some extent, how strongly all particles within are 
correlated with all those outside It is thus, for example, well suited to detect long- 
range correlations near the critical point of a quantum phase transition by enlarging 
see [a El HE]. 

For a many-particle system without long-range interactions the ground-state EE 

H(Q,LVt) = 2S{f),LVt), L>1 (1) 

is widely believed ElE] to grow to leading order proportional to the area \dfl\L^ ^ 
of the boundary surface d{LQ) of the scaled region Lfl as the (dimensionless) scaling 
parameter L tends to inhnity, L —)■ oo. Here, d = 1, 2, 3,... is the spatial dimension of 
H. If the particles are fermions and if there is no spectral gap above their ground-state 
energy in the inhnite-volume limit, the effective long-range correlations lurking in the 
(Pauli-)Fermi-Dirac statistics are expected to slightly enhance the large-scale growth of 
H{0, LQ) by a logarithmic factor ln(L). Indeed, for the free Fermi gas inhnitely extended 
in d-dimensional Euclidean space such a large-scale growth of the ground-state EE 
with a precise and rather explicit prefactor has been proved with full mathematical 
rigour [S], thereby conhrming a stimulating conjecture by Gioev and Klich jH). 

In this note, we present and discuss the hrst rigorous results on the EE of the free 
Fermi gas in in the state of thermal equilibrium at nonzero temperature, T > 0, and 
chemical potential /r G M := The latter we mostly suppress for notational simplicity, 
but also because we often consider thermal quantities for hxed (mean) particle-number 
density p > 0. In contrast with the ground-state or T = 0 case, the EE at T > 0, 
denoted by H{T,Q), must not be expected to be just twice the local thermal entropy 
S{T, G) < oo, since thermal states are mixed ones. However, in trying to dehne H{T, H) 
as the mutual information 




H{T, H) = 5(T, H) + “5(T, W\n)- S{T, M'") 


( 2 ) 


of the bipartition, we are confronted with the two inhnities S{T, M'^\H) = S{T, M'^) = oo 
due to the additivity of macroscopic thermal entropy. We solve this problem by rewriting 
the right-hand side of (|^ in a mathematically and physically reasonable way as a sum 


of two hnite (that is, not inhnite) differences, see (19) and (20) below. By construction, 
H{T, LQ) is then well-dehned and exhibits a leading term proportional to as 

L ^ oo. Given that, our general line of arguments is similar to that of Ref. [TO] 
devoted to noninteracting fermions in the d-dimensional simple cubic lattice with 
emphasis on the case d = 1. 

Our main results may be summarised as follows. For the (spinless) free Fermi gas 
in at any T > 0 we hnd the following two asymptotic large-scale expansions: the 
local thermal entropy satishes 


5(T, in) = s{T)\n\L^ + ri{T, dn)L'^-^ + ... 


(3) 
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H{T,Ln) = 2r]{T,dn)L'^-^ + ... , ( 4 ) 

up to terms growing slower than as L —)• oo. Here, the bounded subregion 

H C may be rather general except that its boundary surface dVL (if d > 2) should be 
sufficiently smooth. For further assumptions see our theorem in Sec. 3 below. There we 
also make the definitions of the entropies S{T, H) and H{T, H) more precise and express 
the coefficient r]{T, dfl) in terms of a multiple integral. Nevertheless, in this Letter we 
concentrate on the physical aspects and publish the somewhat lengthy mathematical 
details in separate papers [IHEIE]. In the next section we just identify s{T) and offer 
some explanations and comments. 

2. Physical meanings of the asymptotic coefficients and their dependence 
on temperature. 

Not surprisingly, the leading asymptotic coefficient s(T)|r2| in ([^ is nothing but the 
thermal entropy contained in hi and s(T) > 0 is the (infinite-volume) thermal entropy 
density or mean entropy. The latter is given by the thermodynamic relation 

s(T) = Tp(T), p(T)~ l^dEM(E)fT(E-iJ.), (5) 

where the integral is the pressure of the free Fermi gas as a function of T (and /i). 
The quantity Af{E) := {27ih)~^ dp Q{E — e{p)) defines the integrated density of 
states AT : M —)■ [0, oo[ of the energy-momentum dispersion relation e : —)■ [0, cxd[ 

which characterises the translation-invariant one-particle Hamiltonian of the free Fermi 
gasjU Here, 2nh is Planck’s constant and 0 is Heaviside’s unit-step function. The 
second factor of the integrand in ([^ involves the Fermi function, /-r : M —)• [0,1], 
/r(F) := [1 -|- exp(F/T)]“^, where from now on we put Boltzmann’s constant ks = 1- 
From ([^ and Sommerfeld’s asymptotic low-temperature expansion, 

ME) = Q{-E) - M/6) TM'{E) + ... (6) 

(in distributional sense), we get the well-known formula 

s(T) = (,TyZ)M'(^)T+... ( 7 ) 

up to terms vanishing faster than T as T j, 0. Eq. ([^ holds, as it stands, at fixed 

chemical potential /i G M. If instead of fi the particle density p > 0 is kept fixed, one 

has to invert the thermodynamic relation p = dp/dfr between p and p,. By ([^, one thus 
finds at low temperatures another well-known formula 

p) = — (tt^/S) -I-... , (8) 

§ For convenience, we have assumed e{p) > 0 for all p so that N{E) = 0 if i? < 0. This is no loss 
of generality as long as e is bounded from below. 
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where Sp '■= limTio ti(T, p) > 0 is the Fermi energy which satisfies p = N'{eF)- 
Consequently, Eq. (0 implies that s(T) = /Z)^'{eF)T + ... for fixed p > 0. For 

the ideal (non-relativistic, free) Fermi gas [H], corresponding to the prime example 
e{p) = p^/(2m) with m > 0 being the mass of each particle, we recall the explicit 
formula 

M{E) = 0(E) [mE/{2nh^)Y^‘^/{d/2)\ (9) 

for its integrated density of states. Quantum effects dominate at low temperatures and 
become weaker at higher temperatures. Accordingly, the properties of the ideal Fermi 
gas approach those of the ideal Maxwell-Boltzmann gas in the high-temperature limit, 
T —)■ cx). For example, in this limit the thermal entropy density of the ideal Fermi gas 
grows, to leading order, proportional to at fixed p and proportional to ln(T) for 
fixed p, in symbols, 

s{T) ~ (T/To)^/' (p fixed), (10a) 

s(T) ~ ln(r/ro) (p fixed). (10b) 

Here, the constant Tq > 0 is an arbitrary comparison temperature. 

Returning to the expansions (|^ and (|^, we note that the other asymptotic 
coefficient, p(T, 911), is also positive. On the one hand, it represents the thermal entropy 
on the boundary surface dLl and determines the first-order finite-size correction to the 
infinite-volume entropy (density). On the other hand, T]{T,dQ)L'^~^ is half the thermal 
EE to leading order in L. Consequently, Eq. ([^ shows that the local entropy at T > 0 
displays a leading large-scale behaviour in agreement with a “volume law” as it should 
be. In contrast, Eqs. (|^ and (|^ show that the EE at T > 0 obeys a strict “area 
law”: although the two subregions Lfl and its complement \ Lfl carry (contrary 
to the case T = 0) extremely different local entropies (namely S{T,LQ) < oo and 
S{T,'R'^ \ LQ) = oo), on their common boundary d{LQ), or interface, the entropies 
are equal and proportional to as L —)■ oo. Roughly phrased, the logarithmic 

enhancement ln(L) present |H1 E] in the large-scale behaviour of the ground-state EE, 
see ([^ and (13) below, disappears when the temperature is raised from T = 0 to 
T > 0 because the Fermi surface “grows soft”. Incidentally, we note that p(T, 914) 
does not depend on the choice of a condition imposed on the domain of the (quantum) 
Hamiltonian at the boundary dQ, because we work from the outset in the infinitely 
extended position space and view all operators to act self-adjointly on the associated 
one-particle Hilbert space L^(M'^) of square-integrable functions "9 : — )■ C, g h->■ ^p{q). 


The coefficient rjfT^dVt) is given by a rather complicated integral (see (22-24) 
below), which, fortunately, is well-known in the theory of semiclassical expansions for 
traces of certain (truncated) Wiener-Hopf type operators, see Refs. [I5l [161 HZl HE]- 
Interestingly, in the limit T 0 the coefficient //(T, dVL) simplifies, displays a logarithmic 
singularity and takes (at fixed /i > 0 [j|]) the rather explicit form 

dO) = (1/12) J(9F^, dO) ln(To/T) + ... (11) 


If /.i < 0, then g{T, dil) vanishes as T | 0. 
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up to terms remaining bounded as T 0. Here, the level set dT^ := {p E : e{p) = fi} 
in momentum space is the (effective) Fermi surface corresponding to p,. The factor 
J(5r^, OLL) is defined as in Ref. [8] and for d > 2 given by the twofold surface integral 

J(9r^,9H) := (27rh)^“'^ f d(j(p)dr(g) |m(p) ■ n{q) \ . (12) 

Jdr^xdQ 

The vectors m{p),n{q) G denote the exterior unit normals at the points p G c^T^ and 
q G dfl, respectively. The canonical {d — l)-dimensional area measures on the surfaces 
dT^ and dfl are denoted by a and r, respectively. If we fix the particle density p > 0 
and use ([^ in 0 , we arrive at 0 with p replaced hj ep- 

By identifying the large ratio Tq/T inside the logarithm in ( pT| with the scaling 
parameter L, Eqs. (|^ and ([^ give 

S{0, LQ) = (1/12) J{dr^, dQ) ln(L) + ... (13) 


in agreement with the result for T = 0 in Refs. |Hl[9] (resp. the corresponding expression 
with p replaced by ep). For an isotropic dispersion relation e we know from Ref. [S] 
that J(cIF^, dVt) is proportional to the area This is even true for p(T, dO) itself, 

at arbitrary T > 0. However, the emerging prefactor, the thermal entropy surface 
density, remains to be given by a multiple integral, see the remarks below the subsequent 


Eq. (24) 


As for the entropy density s{T), the leading high-temperature behaviour of the 
coefficient p(r, dLl) of the ideal Fermi gas depends on whether /i or p is kept fixed. 


More precisely, it follows from (24) that 


7 ]{ T , dn ) ~ ( t / To )(‘^-^)/2 
p(T,aH) ~ (To/T)1/2 


(p fixed), 
(p fixed). 


(14a) 

(14b) 


as T —)■ oo. Eq. (14b) reflects the fact that the particles become effectively uncorrelated 


for fixed particle density at sufficiently high temperature. 


3. Precise definitions and formulations of results. 

In order to define the local thermal entropy and the thermal EE precisely we first recall 
that the infinite-volume equilibrium state of the free Fermi gas at temperature T > 0 
and chemical potential p G M is quasi-Gaussian (in other words, quasi-free) and uniquely 
determined by its reduced one-particle density operator fT{e{P) — pi) on L^(M'^). Here, 
P := —ihd/dq is the canonical-momentum operator, e{P) > 0 the one-particle quantum 
Hamiltonian and 1 the identity operator. The local (or truncated) version 

D{fp,n):=lnfT{e{P)-pl)ln (15) 

^ In particular, for the ideal Fermi gas in one simply has J(dr^,9n) = 2A/d-i(/i)|dB|, where 
Nd-i{E) is given by the right-hand side of § with d replaced by d — 1. 
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of the density operator then characterises the qnasi-Gaussian substate obtained from 
the equilibrium state by spatial reduction to G C Here, —)■ 

denotes the projection operator associated with the indicator function of that is, 
{lQ^p){q) := fj{q) if g G G and 0 otherwise for all E L^(]R'^). Next we recall the binary 
entropy function h : [0,1] —)■ [0, ln(2)] dehned by 

h(0) := h(l) := 0 and h(t) := —t In (t) — (1 — t) In (1 — t) if fG]0,l[. (16) 


The positivity 

h{X,r,t) := h((l — A)r + At) — (1 — A)h(r) — Xh{t) > 0 (17) 

for all X,r,t E [0,1] is equivalent to the concavity of h. In fact, h is even operator- 
concave [19]. The local thermal (von Neumann) entropy is now given as the trace 


5(T, H) := Tr[ln h(Zl(/T, G))lo] = Tr[ h(Zl(/T, G))], 


(18) 


see, for example. 


The second equality in (18) follows from h(0) = 0. 


The dehnition of the thermal (von Neumann) EE requires two steps. In the hrst 
step we introduce the “entropic” operator difference 


A(T, a) -- lnh{D{fT,a))l» - D(h o h,a). 


(19) 


Here, the operator D{h o /t,!!) is obtained from (15) by replacing the Fermi function 
fx with the composed function ho fr dehned by {h o fT){E) := h{fT{E)) for all E eM.. 
From the operator concavity of the function h and from Refs. |T9[|2T] we get the operator 
positivity A{T,Q) > 0. In the second step we dehne the thermal EE for a bounded 
H C as the sum of two positive traces 


H{T, n) := Tr A(r, G) + Tr A(T, W^\n). (20) 

This is the precise version of (|^. Arguments as in Ref. [22| show that even the 
second trace is hnite, although the (positive) minuend and the (positive) subtrahend of 
A(T, \ H) have both an inhnite trace. Therefore, we arrive at the (in)equalities 

0 < H{T,n) = H{T,W^\Ll) < oo. ( 21 ) 


Moreover, in the limit T j, 0 we get back to Eq. ([^ by observing that h o /q = 0, confer 

(§. 

To explain the coefficient r]{T, dQ) of the subleading asymptotic behaviour of 
S{T, LQ) and of the leading behaviour of H{T, LQ) as L —)• oo, we need some auxiliary 
dehnitions. In contrast to the well-known leading “volume term” in ([^ the subleading 
“area term” is new and rather complicated. It cannot be obtained from simple heuristic 
considerations, not even for d = 1. But it can be derived from the semiclassical “area 
coefficient” in Refs. ng usi [m ng by observing that the traces in ( jig ) and ( [20| , 
with replaced by LQ, depend on Planck’s constant 2nh and the scaling parameter 
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L only via the ratio h/L which can be seen by a (unitary) dilatation. In order to 
recall this “area coefficient” from the mentioned literature we first define the function 
U : [0,1] X [0,1] [0, cx)[, (r, t) ^ U (r, t) by 




h{X,r,t) 


see (16) and 0 for the definition of h. Then we consider the integral 

U{g{v),g{w)) 


U\q\ := lim 

^ ^ <5i0 


dn dw 0(|n — w| — 5) 


( 22 ) 


(23) 


{v — wY 

defined, in the principal-value sense, for smooth functions : M —)■ [0,1]. We observe 
that hl[g] >0 since U{r, t) > 0 due to the concavity of h. If the function h were smooth, 
then U[g] could be defined as a standard Riemannian integral, whose finiteness could be 


easily checked (for smooth g). For the function (16) however, as well as for certain other 


functions (such as f h->• f“(l — f)", a G ]0,1[ ) being continuous but not differentiable at 
the points t = 0 and t = 1, the finiteness of L([g] is a non-trivial matter. This and other 


relevant properties of the integral (23) are investigated in m- 

If d > 2, then we consider for g at given T > 0 the 2{d — l)-parameter 
family of functions fT-,{p,q) : R —t [0,1] defined in terms of the Fermi function by 
fT-,{p,q){v) ■= /t(^(p + vn{q)) — /i), where n G M, g G 911 and p G T*(911) with 
T*(9r2) = being the dual space of the {d — l)-dimensional tangent space of dVt at 
the point q. The vector n{q) G and the subsequent area measure r have the same 


meanings as in (12). Finally, we define 


r]{T, dfd) := (27rh) 


l-d 


/ dr(g) / dpU[fT-,ip,q)]. 
'an JT*{an) 


(24) 


If the dispersion relation e is isotropic, then the functions fT;{p,q) do not depend on the 
parameter q G dVt due to the orthogonality of p G T*(9r2) and n{q). Consequently, the 
surface area |9r2| can be factored out on the right-hand side of ( |24| ). Nevertheless, the 
multiple integral underlying ri{T, dfl) remains to be a fourfold one for d > 2. For d = 1 
the set n C M is the union of finitely many pairwise disjoint bounded intervals, according 
to an assumption in the subsequent theorem. Then |9r2| equals the (even) number of 
all endpoints of the constituent intervals and one has r](T,dfl) = U[fT o (e — p)] |9r2|, 
which still involves a threefold integral. 

Now we are prepared to state our central “technical” result as the following 

Theorem. For any temperature T > 0, for any bounded subregion hi C with finitely 
many connected components and (if d> 2) piecewise smooth boundary surface dfl (see 
Ref. fMBf) and for any smooth and polynomially bounded dispersion relation e the trace 
of A(T, A) is finite and positive for A = r2 or A = M'^\r2 with the same leading term in 
both large-scale expansions 

Tr A(T, LA) = r/(T, dLl) L^-^ + ... , 


up to terms growing slower than L^ ^ as L —?• oo. 


( 25 ) 
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Our proof goes as follows. At first we “smooth out” the function h, which enables 

or |T7] 


us to refer for the formula (25) directly to available literature, notably to 


(see, in particular. Chap. I and V in irzi)- To return to the original non-smooth function 
h, we then “close” the asymptotics using the estimates obtained in [22] for non-smooth 
functions of operators of the type D{f,Q). This strategy is similar to the one applied 
. Full mathematical details will be published elsewhere m [T3|. There we also 


m 


establish the validity of formula (25) even if L —)■ cx) and T 0 simultaneously provided 


that LT/Tq > 1 throughout. In our view, these observations provide a deeper insight 
into the low-temperature behaviour of the thermal EE. 

Combining (20) and (25) immediately gives the claimed large-scale behaviour (|^ 


of the thermal EE. From (25) we also infer the claimed two-term large-scale behaviour 


(|^ of the local thermal entropy (18) by observing 

5(T,Lf2) = TrT)(ho/r,Lfi) + TrA(T,Lfi) 
= s(T)|f2|L'^ + TrA(T, LLl) 

= s{T) \Q\L^ + 7]{T, dQ) + .. 

where we have used the identities 


{2iih) 


-d 


<iph(Me(p)-li))= / dE^f'(E)h{fT{E-f,))^s(T) 


(26) 


(27) 


in the second equality. The leading large-scale behaviour of the local entropy S{T, LQ) 


at T > 0 was hrst proved in 


The subleading correction of the order L 


d-l 


m 


(26) is new. 


4. Summary, discussion and an open problem. 

For the free Fermi gas in multidimensional continuous space in thermal equilibrium at 
temperature T > 0 and for a given bounded subregion hi C we carefully distinguish 
between the local thermal entropy S{T,Q) > s(T)|r2| and the thermal spatially bipartite 
entanglement entropy > 0. Their large-scale behaviours ([^ and (|^ contain the 

thermal entropy density s(T) > 0 and the new asymptotic coefficient r]{T, dQ) > 0 as 
two characteristics of the rather old free Fermi-gas model [Tl|. The results (|^ and (|^ 


are physically even more relevant than the corresponding ground-state results (13) and 


([^, because in real gases the temperature is never strictly zero. Furthermore, the new 
results ([^ and (|^ deepen our understanding of the older ones (13) and ([^ by observing 
that s(0) =0 and that t]{T, dQ) diverges logarithmically as T 0 according to ( 0 . A 


result similar to (11) for noninteracting fermions in the one-dimensional lattice was 


derived in ra. but without an explicit prefactor. As in [8], HOj, many of our present 
results extend to the whole one-parameter family of (quantum) Renyi entropies, see 

[niiis]. 

Finally, we check whether the exact asymptotic results (13), 0, 0 and (0 for 
the local thermal entropy S(T, LQ) can be descibed consistently by a so-called universal 
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crossover formula as obtained in [26l EZl EHl EHl EQI [3T] by arguments from a suitable 
conformal field theory, at least for d = 1. Such a formula has the appealing form 


S(T,La) = L‘‘-Mln 


— sinh 

/ T\1 

L—] 

Lt 



( 28 ) 


with suitable constants A > 0 and Tq > 0 not depending on L and T. Remarkably, if 
we choose A = (1/12) J(dr^,dQ) and Tq = ?)A/[n‘^N''{pi)\Vt\\^ then we hnd consistency, 
provided that L S> 1 and T <^Tq. We conclude that ( |2^ correctly reflects asymptotic 
properties of the free Fermi gas in if L is large and T is small, but does, for example. 


not reproduce the large-T behaviours (10) and (14) (of the ideal Fermi gas). However, 


Eq. (|28|) suggests the scaling-type formula 

- Sifd.LVt) 


li™ TXT 

L^oo ^ 


5 


(xTq 


I 


= H In 


sinh(a;) 


X 


( 29 ) 


for any ;r > 0. Although we have a result for the simultaneous limits L —)■ cx) and T 0 
with X = LT /To (> 1) kept hxed, see a remark below the above theorem, at present we 
do not know whether (p^ or a similar formula follows from the “microscopic” dehnition 


(18), not even for d = l. 
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